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Abstract 

We discuss the problem of CMB spectrum corruption during Galac- 
tic emission removing. A new technique of spatial-frequency data re- 
duction is proposed. The technique gives us a possibility to avoid a 
spatial harmonics nonorthogonality. The proposed technique is ap- 
plied to the two-year COBE DMR sky maps. 

We exclude the harmonics with 1=7, 9,13, 23 and 25 as having 
anomalous statistics noise behavior. One shows that procedure do 
not give systematic errors, if the data are statistically regular. 

The spectral parameter of the power spectrum of primordial per- 
turbation n = 1.84 lb 0.29 and quadrupole moment Q2 = 15.22 it 3.0 
are estimated. The power spectrum estimation results are inconsistent 
with the Harrison-Zel'dovich n = 1 model with the confidence 99%. 
It is shown a necessity of an increasing a survey sensitivity to reach a 
more reliable estimation of the cosmological signal. 



*accepted to Pis'ma v Astronomicheskij Zhurnal, 1996, v. 22, in Russian 







1 Introduction. 



A large scale cosmic microwave background (CMB) anisotropy carries the in- 
formation about a primordial metric perturbation amplitude and spectrum, 
it shape reflects the Universe evolution and structure. To encrypt the in- 
formation one analyzes a CMB anisotropy spatial spectrum and compares 
it with a theoretical models with different types of the Universe evolution 
scenario. A sensitivity rising in the CMB investigations caused a detection 
of the anisotropy. As a result one may find prerequisites of a more detail 
examination of the physical conditions involved in the Universe birth and 
evolve. Being high sensitive and multi-frequent COBE experiment (Smoot 
et al., 1992, Bennett et al., 1994) gives a possibility to careful expanding of 
our knowledge about the CMB structure and spectrum. 

Unfortunately, the Galaxy radiation predominates even at a millimeter 
wavelength range. This rises a problem of the radiation filtering. The 
most Galaxy radiation is coming from rather narrow region located along 
the Galactic plane. The main method to exclude the radiation is to "cut 
out" the Galactic plane vicinity and to analyze the rest regions on the sky 
map ( Strukov et al., 1991a, 1991b, Smoot et al., 1992, Bennett et al., 1994, 
Gorski, 1994, Gorski et al, 1994, Wright et al., 1994b). This method is well 
known but it works perfectly only if variance data analysis is applied. If 
we use the Galactic "cut out" in the case of spatial spectrum analysis we 
lose a spherical functions orthogonality on the rest sphere part. As a result 
the harmonics begin to influence each to other. The problem may be solved 
partly using a new basis which is orthogonal over the rest part of the sphere 
( Gorski, 1994) or is orthogonal only to a monopole and dipole because the 
excluding of these two components effects the most spectrum modification. 
The spectra obtained this ways may differ from a real one (Bunn et al., 1994) 
and may cause the essential different conclusions when interpreting the re- 
sults. 

A modern multi- frequency experiments seemingly allow us to separate the 
blackbody CMB radiation component from the frequency dependent Galaxy 
foreground. But variations of the radiation spectral composition prevent the 
procedure from being performed with a desired accuracy. 

We propose the method of the spatial spectrum analysis which allows to 
reach the most spatial resolution without any loss in the reciprocal harmonics 
orthogonality. The method imphes that the CMB filtering stage is conducted 



1 



not before but after the radio map is expanded to spherical harmonics. 

There are two filtering stages in the proposed process. The first stage is a 
spatial filtering one when some components related with the Galactic plane 
and the Galaxy center are excluded from the total spatial spectrum. By this 
expedient the Galactic radiation magnitude is essentially reduced. And only 
after that, during the second filtering stage a "frequency cleaning" procedure 
is used. At this stage the CMB anisotropy is separated from the Galactic 
radiation using several frequency channels and some a priori Galaxy radiation 
model. The both stages nowhere are associated with any signal harmonics 
orthogonality losses. 

The main method difficulty may be found at the second stage and is 
bound up with the problem of a correct Galactic radiation model selection. 
But using the first stage filtering we strongly reduce the effect of the model 
uncertainty on the net result. 

As an initial data we used public released COBE DMR 2-years all sky 
maps from 31.5, 53, and 90 GHz frequency channels. In addition we exploited 
a 19.2 GHz balloon data kindly provided by the authors of the survey. 

2 The Method. 

A data reduction process practically always is accomplished by some types of 
the data censoring like weighing, excluding outliers, corrections for systematic 
errors, trends and so on. A censor method selection is based on a priori 
information about the analyzing sample. The censoring is aimed at the 
improvement of the estimation stability and at the approaching the real data 
statistic parameters to the used a priori models. 

In our case the values on the observed radio map is a sum of a cosmo- 
logical signal we try to measure, a receiver noise, the Galactic radiation, and 
perhaps some signal related with neglected systematic effects. It is the cos- 
mological signal that we are interested in, so the other components remain to 
be carefully estimated, significantly reduced or to be completely excluded. In 
an ideal case (no Galactic radiation and no systematic effects) the procedure 
would not cause a significant change in the cosmological signal parameters. 

In a Galactic coordinate system one may find a Galactic emission concen- 
trated near the Galactic plane and around the Galaxy center. The radiation 
may be described using only a few spatial harmonics. If these harmonics 
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are excluded the Galactic radiation is substantially reduced. In our case we 
excluded 23 harmonics only. At the same time the rest spatial spectrum 
consist of 649 harmonics. If a signal is normal distributed the 649 harmonics 
are quite sufficient for a stable signal spectrum parameter estimation. 

To estimate the signal parameter adequately it is absolutely necessary to 
know (or to estimate) exactly a noise which may be found in any analyzed 
data. The most problem in the procedure is the presence of anomalous mag- 
nitudes in the noise spectrum. If the signal/noise ratio is not good enough 
the anomalies may cause a significant variation in the parameters to be esti- 
mated. 

We emphasize that a cause of the anomalies may be some residual sys- 
tematic errors. In this case it is very difficult to build up a statistical model 
of the anomaly and the accuracy of the signal estimation may drops below a 
point that can be tolerated. 

If the number of the anomalous harmonics is not large in comparison 
with the total signal harmonics, the radical solution of the problem is an 
exclusion of the anomalous harmonics from the analysis. If data are regular, 
the excluding some anomalous (in the noise spectrum) harmonics is rather 
safe procedure. That is because the cosmological signal and the noise are 
Gauss type, and the signal parameters are weakly depended on both the fact 
of excluding and the criteria of anomalous outliers selection (in our case we 
rated to anomalous the outliers beyond the 3-sigma). 

But if the neglected systematic errors are present in the data, this proce- 
dure causes the significant improvement of the sample regularity and causes 
the increasing of an estimation efficiency. At the same time if the excluding 
process significantly changes the estimated parameters, it may be a strong 
reason to assume the data is not regular, and consequently may indicates 
that it is necessary to use some type of procedure to reduce the irregularity. 

It should be noted that using data reduction methods which break a 
spatial harmonics orthogonality causes a certain spectral smoothing and may 
reduce or even may camoufiage the abnormal amplitudes. In this situation 
the risk of missing a really existing anomaly is run, and the signal estimation 
process may be accomplished with some additional errors. 

The data reduction procedure we propose is the following. We use COBE 
A and B channels (Wright et al., 1994a) and produce the sum and the dif- 
ference radio-maps as a half-sum and a half- difference of the A and B maps 
respectively. All data on the sum map are a sum of a signal and a noise. All 
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data on the difference map are only the noise, the signal amplitude is zero. 

We exclude from the sum map the monopole and dipole components, and 
the spectral components cased by the Galaxy emission. Then the sum map 
is additionally "cleaned" from the Galactic radiation. For this purpose we 
use 19 GHz and 31 GHz maps to determine the Galactic radiation spectral 
index. Then the 31 GHz map is subtracted from 53 GHz map, but with a 
weight corresponding the determined spectral index. 

We assume the hypothesis the cosmological signal and the noise are nor- 
mal distributed. The noise spectrum shape is assumed to be completely de- 
termined by the number of measurements in every point on the sky. Based on 
this we calculate the noise parameters and discover anomalous noise compo- 
nents. Then we use a Monte Carlo simulation to estimate signal magnitude, 
signal spectrum parameters, and the estimates confidence levels. The latter 
procedure is applied both to a total analyzed spectrum and to a censored 
spectrum with anomalous components excluded. In addition we test the de- 
termined estimations for stability to the procedure of quadrupole component 
subtracting. 

2.1 Signal Analysis on the Sphere. 

The signal on the sphere may be represented as a spherical harmonics ex- 
pansion Yi^rn{0,'P)- 

AT{0,ip) = ^aTYUO,v) 
In this case the signal variance on the sphere is represented as: 

^'^ I m I 

where AT^ is a /-th spherical harmonic power. 

In any real experiment the analyzed signal af^signai affected by a transfer 
function Wi (thereafter we use the COBE antenna transfer function (Wright 
et al., 1994a) as Wi) and some noise a^^^^ise is added: 

AT{e, <p) = Y.^Wla^s^,nal + «rnoi J>^/,-(^. 

The cosmological signal on the sphere is a sample of a random process. 
In a standard cosmological scenario the inflation creates adiabatic scalar and 
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tensor Gauss fields. So the spectrum is assumed to be a sample of a normal 
random process with independent harmonics, (a™a^) ~ 6{l, k)6{m,n). 

If initial density perturbation spectrum is {5p/ pY = Ak", one may find 
for the mean lowest multipoles in expansion AT{6, (f) (Bond, Efstathiou, 
1987): 



r(/ + (n-l)/2)r((9-n)/2) 



(«) ) ~ ^ (5 _ n)/2)r((3 + n)/2) 

A consideration of other physical phenomena ( evolution effects, acoustic 
oscillation. Silk effect etc.) causes a more complicated spectrum (White et 
al., 1994 and references therein). Hereafter we approximate a real spectrum 
with the (H). For a quantitative spectrum description we use only two para- 
meters: a power spectrum index n and quadrupole magnitude Q^ms — {^'^2) 
as an amplitude. 

A correlation function on the sphere C(/5) = ((AT(gi))(AT(g2))), ^1^2 = 
cos(/3) may be represented in terms of spherical functions: 

^(/^) = ^E^Kcos /3)E(an^ 

Using the spherical harmonic expansion we may calculate the correlation 
between two maps (for example between 31.5 GHz and 53 GHz maps): 

p(MAPl,MAP2) Ei,rnar{l)ar{2) 



v/D,™(ar(l))'E/,„.(ar(2))2 
2.2 Spectrum Parameters Analysis. 

A cosmological model scenario gives us a prediction of a mean CMB spectrum 
fluctuation. Experimenters deal with the single sample of the scenario. In 
addition a noise is always present in a real data. In this connection there is 
the problem of a signal parameters estimation. The problem was discussed 
by Bunn et al., 1994 and Sazhin et al., 1995. Here we will not go into the 
problem, but will touch upon the estimation not only the signal, but the 
noise parameters too. 

The maximum likelihood estimation is known to be the most effective. It 
may be illustrated if the signal spectrum has a given a priori shape. 
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Let the spectrum shape is ((a™)^) = Q'^Fim- If the amphtude is Q, the 
hkehhood function may be described as: 



The maximum hkehhood estimation Qml may be found if solve the equa- 
tion maximizing the hkehhood function: 



dlnfiQlFi] 



2 



dQ 



l,m 

where M is the number of the analyzed harmonics. 

The parameter is evident having Xm distribution with M degree of free- 
dom, Qml ~ Xm- The amplitude estimation calculated using the variance 
on the sphere is: 







QpowER - j2{AT^)\Q = l ~ al\Q = 1 ^ 



It is clear that: N^ff = . Jas. < M. So for the Harrison-Zel'dovich 

spectrum iVg// ~ 100, M ^ 700 p Imax = 25. The maximum likelihood 
estimation accuracy is about 5 % and the total power estimation is about 
14%. 

Using the approach we may determine the signal spatial spectrum para- 
meters in the presence of noise. Maximum likelihood estimation offers to find 
a minimum of the functional: 



relative Q and n with given spectrum shape and measured signal and noise 
sample spectra. Here {{ai{Q,n)Y) is the model signal spectrum, depending 
upon desired parameters and described as (0), (&[")^ is the measured spec- 
trum (we describe it further), and {{ciY^noiseY) is a model noise spectrum (|^), 
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its amplitude is found using the same approach but from the difference map 



In addition we are to take into account a procedure of the Galactic radi- 
ation "cleaning" and a procedure of a conversion the antenna temperatures 
into a thermodynamic scale. In all cases the summation is taken over mul- 
tipolcs chosen for the analysis. The necessary conditions for the described 
procedure is a knowledge of the measured signal and noise spectra. 



There are a lot of papers concerning the COBE data noise analysis. Bennett 
et al., 1994, Lineweaver et al., 1994 discussed the correlation noise properties 
and showed a weak ( 0.45% ) correlation at 60°. So in practice we may 
assume the noise is uncorrelated. The noise power estimation may be done 
using the difference maps. 

As a first approximation assume that the noise is white. However the 
real observations have different accuracy on the sphere. The COBE orbit 
configuration was so that the points near the ecliptic poles were observed 
much longer than others. Let us consider the methods to take into account 
the measurement accuracy is not equal for different points. 

Let the noise in every point on the sphere be normal with zero mean 
and the variance a^{9, (p) = a^/N^O, (p), where N{9, ip) is a number of mea- 
surements in the point with the 9, (p coordinates. Then we may propose 
the following noise spectrum model on the sphere. If there is a white noise 
sample G = ^n^Yl^ ^ {ji^vJ^ ) = 0, {{n^Y) = n'^ = const, the true noise on 
the sphere is a product of G{9, ip) (here, and only here n is a magnitude of 



the noise power spectrum, not a spectrum index), and 1/ JN{9,ip): 



2R{9, if) = l/^N{9,ip)G{9, ip) = w{9, ip)G{9, ip) 
According to Peebls (1980), assume w(n)Y/^{n) = Ei'm' wJ^'^'Yv"' {^), 



2.3 COBE Data Noise Analysis. 




where 
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Taking into account noncorrelatedness of white noise harmonics a cova- 
riation of a[^oj^g may be written as: 

Wl^noise^^V', noise) = 51 X] '^w"^ ^P/""* {^'v' ^T' ) 
I'm' l"m" 

And finally noise power spectrum is: 

iKno^sel') = n'Y. = I dnw\n)\Yrm' (3) 

I'm' •' 

We use not only unit or zero weights, but unrestricted values. Because of 
that our equation differs from Peebls's one in the squared weights. 

Then, using maximum likelihood approach to the difference maps we may 
calculate the amplitude of the spectrum given above. 

The power spectrum model make it possible to analyze the noise spectrum 
measured from the difference maps. The analysis shows that in general the 
noise spectrum is well described both the model @ and the white noise 
model. It is due to a rather perfect COBE coverage of the celestial sphere. 
But the same analysis shows for the both models the existence of harmonics 
with an excess power. So in the noise spectrum of 31 GHz map one may find 
for / = 13 the power is greater than predicted from (j^). The power excess is 
3.4cr/=i3, where predicted multipole power variation ai is calculated basing on 
the model and on the degrees of freedom for the analyzed component. Further 
we use a censoring of the data to exclude this and similar components from 
the analysis. 

The noise spectral component with / = 25 has anomalous high correlation 
between 31 and 53 GHz maps. After subtracting the Galaxy radiation the 
excess power is 3.3(7^=25. Components with / = 25 were excluded from the 
analysis too. 

We find a significant power excess of several noise components on 53 
GHz map. We exclude the components with power excess more than 30"^. 
After that we determine the spectrum parameters for the rest components 
and the procedure is repeated until the spectrum is free from the abnormal 
components. Thus we excluded the components with / = 7, 23, 9 having the 
initial power excess of 3.11(Ti=7, 2.63(Tz=23, and 2.31cr/=9. 

In the final analysis we use spectral components from l^in = 2 to l^ax = 
22, except components / = 7, 9, 13. 



I'm' 
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Such censoring effects only a weak power spectrum distortion if the signal 
amplitude distribution is normal. The matter is discovered by numerical 
modehng, the results are presented in Chapter 3. 

During signal spectrum analysis we ignore a noise spectrum components 
correlation and use the spectrum shape only. Further modeling shows the 
same results if we use a pure white noise. 



2.4 Selection of the Galactic Radiation Model 



Our approach is based on the spatial data filtering which reduces the in- 
fluence of the Galactic spectrum index uncertainty. The accepted Galactic 
radiation model for the analyzed frequency region is assumed to involve a syn- 
chrotron component with frequency dependence as T = Tqu", a = —3 ± 0.2 
and a bremsstrahlung component with spectral index a = —2.1 (Bennett et 
al., 1992). If there are the space regions of a different nature of radiation on 
the line of sight, it may cause a spatial variation of the spectral index a. Let 
us try to estimate roughly the contribution of such variation if we use some 
effective but constant value ao. Let the true value of the spectral index be 
a = ao + Aa. After the frequency depended part is excluded from the data 
Ti T2 of two frequency channels, one may derive: 



6T 



T, 



1^1 



T2 Aa In ( — 



The residual signal power caused by the spectral index variation ay^j^ = 
AQ;^(ln z/2/z/i)^cr^^P2. For the frequencies 31.5 and 53 GHz the residual RMS 
is smaller than 

0.52Aa(TMAP2- It can be seen that the less is Galactic radiation contribu- 
tion, the less is the residual RMS on the more high frequency radio map. 

It is evident the most Galactic radiation is concentrated near the Galaxy 
plane and near its center. If the spatial spectrum is represented in Galactic 
coordinate system the most power may be found in the components with 
m = for the even multipoles (Galaxy plane) and with m = 1 for the odd 
multipoles (Galaxy center). So if we exclude this components the influence 
of the spectral index variation may be reduced drastically. At the same time 
this procedure does not break the multipoles orthogonality. 

Table |^ shows some results obtained after the spectrum analysis from 



^min 2 to 



25. 
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In addition one may find at the last row the difference between 53 GHz 
map and the 31 GHz map scaled with a spectral index a = —2.15 after 
spatial filtering. Data are represented in antenna temperatures. One can see 
the resultant asky niay be explained either by a cosmological signal existence 
or by the spectral index variation within Aa ~ 0.5. 

Unfortunately the available data give no way to estimate the variation 
with required precision. For more reliable conclusions it is necessary to have 
the sensitivity many times better than COBE instrument has. 

We tried to estimate the spectral index variation crude. For this purpose 
we use both the COBE data and 19.2 GHz survey (Boughn et al., 1992). We 
has converted latter data to a COBE beam shape (Wright et al., 1994a). Af- 
ter that we analyzed the frequency dependence of the most intensive spherical 
harmonics and determined the accuracy of the estimation. 

The data are shown in Table |I] where one can see correspondingly: ana- 
lyzed frequency ranges, spectral index a for desired signal/noise ratio equal 
to 5, spectral index variation Aas, and the predicted spectral index variance 
AaN, calculated basing on the noise analysis. In addition there are shown 
the spectral index aw and its variation weighted-mean for all harmonics with 
m 7^ for Z = 2A; and with m ^ 1 ioi I = 2k + 1. 

It may be noticed the absolute value of the spectral index ai9_3i is regu- 
lar lower than a3i_53. The contribution to the total spectrum of the spectral 
components with signal/noise ratio equal to 5 is more than 40/found the 
only spectral component with the spectral index like a synchrotron Galac- 
tic radiation has. This is 02,2 component with the spectral index aig-si = 
-2.86 ± 0.15, a3i_53 = -2.92 ± 0.22. 

The different maps correlation analysis gives the mean spectral index 
estimate in the region a = —(2.2 2.3). 

So we may assume the hypothesis the spectral index is constant and its 
measured variations are caused only by the instrumental noise. 

On the basis of the analysis we can conclude for the most significant 
Galactic spectral components (with the exception of 02,2 ) the Galactic radi- 
ation model may be assumed to be one-component, with the radiation cased 
only by an ionized hydrogen. So the Aa variations may be taken as zero. 

Now we may determine the measured spectrum 
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br = (m53) - k{31) (^) (ar(53 GHz) - ar(31.5 GHz) (^) 

where /c(31) = 1. /1. 025724, A;(53) = 1./1.074197 are the scale coefficients 
to transform the thermodynamic temperatures to antenna temperatures for 
31.5 GHz and 53 GHz, a]^{31.5GHz), aY'idSGHz) are multipole coefficients 
for 31.5 GHz and 53 GHz maps. 

Unfortunately, the COBE sensitivity is not good enough and the accuracy 
of the estimation is determined by the instrumental noise. We are pinning our 
hopes on the "Rehct-2" space experiment. "Rehct-2" instrument is preparing 
now, it has an order better sensitivity than COBE has, and it will give the 
possibility to analyze the Galactic spectral index radiation with high relia- 
bility and so to separate the cosmological signal from the Galactic radiation 
with high accuracy. 

3 Results 

As a result of the analysis we obtained the spatial spectrum parameters 
estimates under the assumption that Galactic radiation is one-component. 
For the quadrupole we found: 

Q2 = 15.22 ± 3.0//K 

and for the power spectrum index: 

n = 1.84 ±0.29 

We analyzed spectral components from 1^^^ = 2 to l^ax = 22 if m 7^ for 
I — 2k, and m 7^ 1 for Z = 2A; -|- 1, and (due to the noise anomahes) if 
I ^ 7,9,13. 

The accuracy of the obtained parameters is derived using Monte-Carlo 
simulation. For all estimates we show RMS errors. The procedures of dif- 
ference and sum maps producing and of noise parameters determining are 
included in the Monte-Carlo simulation. The simulation also takes into ac- 
count the noise, its variations, and a cosmic variance. 
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We find the total power measured in chosen spatial frequency window as 
(70.28 /iK)^, noise estimate as (56.76 /iK)^, and on-sphere signal estimate as 
(41.44 /iK)^ ± (14.68 /iK)^. Total number of analyzed harmonics is M = 446. 

Table ^ shows the results. In addition we include into the Table ^ the 
data obtained with anomalous harmonics, and with excluded quadrupole. 

We tested an effect of the harmonics excluding procedure on the signal 
estimation. In the case of the noise outliers caused by systematic errors and 
if all harmonics are orthogonal the excluding does not effect the estimation 
bias, but may only increase the errors in comparison to non-censored data. 
Simulation shows the Q2 and n variations are 0.9 and 0.11 respectively. 

On the other hand, if the outliers are pure stochastic (i.e. we are dealing 
with a low-probability noise sample) the noise estimation bias may arise and, 
as a result may arise the signal estimation bias. We simulated this occasion 
for the signal and noise with the given spectrum and to examine noise outliers 
stronger than 2.3(j. After that we calculated the signal estimates both for 
censored and non-censored spectra. 

After 1228 simulations we have found 500 events with the outliers. The 
corresponding estimates are: {ucensor - n juii) = 0.07±0.08 {Qcensor-Q fuii) = 
—0.4 ± 0.83. It should be noted that the obtained estimates really are upper 
bounds rather than two-sided limits. It is because we exclude the compo- 
nents, which then are used in some linear combination and thus are addition- 
ally normalized. It may, in general, reduce the effect of spectrum censoring. 

COBE data analysis shows (see Table that the procedure of the anoma- 
lous harmonics excluding effects the significant change in n. If the instrumen- 
tal noise is completely random the probability to obtain so large difference is 
smaller than 0.02%. One may see that the applying the anomalous harmon- 
ics censoring to COBE data causes a dramatically decreasing of the strange 
dependence of the results on whether or not the quadrupole is excluded. It 
may be an indirect support to the necessity of proposed data censoring. 

Basing on the obtained parameters we may predict the signal in an ex- 
periment like it is conducted in Tenerife (Hancock et al., 1994). Assuming 
5.5° beam, 8.1° antennae separation, and 3-point method of observation i.e. 
AT = To — 0.5(Ti + T2), we may calculate the following data for the spec- 
trum ([H) type using previously obtained Q2 and n, excluding components 
with anomalous noise and applying the spatial filtering: aTenerif = 54.82 /xK. 

The same analysis, but without excluding anomalous components shows 
two estimates: 
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CTenerif = 37.22 /iK with quadiupole included, and: (JTenerif = 34.47 jjK 
with quadrupole excluded. 

We do not use a 4° binning in the analysis. The binning may decrease the 
obtained data in some degree. The binning is used in Tenerife experiment 
and the results are (Hancock S. et al., 1994): 

(Ts = 49 ± 10 fiK for 33 GHz channel, and = 42 ± 9 fiK for a sum of 
the 15 and 31 GHz channels. So the cosmological signal spectrum with the 
parameters we obtained is in a good agreement with the observations which 
are more sensitive than COBE for high spatial harmonics. 

4 Conclusions. 

The systematic effects exclusion problem is met practically in any experi- 
ment. In the case of CMB anisotropy observation the Galaxy radiation is 
the primary effect. The method usually used to suppress the radiation is to 
"cut off" the Galaxy plane region. It causes the problem of the monopole 
and dipole exclusion. The latter in its turn effects the additional systematic 
errors but at the less level. 

The orthogonal basis when coupled with the spatial-frequency filtering 
of the Galactic radiation allows to avoid the systematic errors mentioned 
above. In addition the proposed approach make it possible to analyze and 
to exclude more fine effects usually caused errors in spectrum parameters 
measurements. If an instrumental noise is normal the data censoring could 
not cause a significant difference between the spectral parameters derived 
from censored and non-censored data sets. Being detected the difference 
may be result from either neglected residual effects or a noise nonnormality. 
Being in the context of the normal distribution function we are forced to 
exclude the spectral components with an anomalous noise behavior. 

The cosmological spectrum parameters we obtained differ in some degree 
from the parameters derived in previous works used the same initial COBE 
data. The difference is due to the anomalous noise harmonics exclusion rather 
than spatial-frequency Galactic radiation filtering. In its turn the harmonics 
exclusion is possible because of the harmonics orthogonality conservation. 

A number of investigators (Table ^ announced the strange end result 
sensitivity to a quadrupole component. Most likely it is also attributed to 
anomalous noise harmonics influence. 
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The results based on 2-year COBE data and obtained by several authors 
are shown in Table ^. In addition we would like to remind the result of 
n = 1.7, derived by (Hancock S. et al., 1994) after a comparison of 1-year 
COBE and Tenerife data. 

So we obtained the estimate of n more precise and somewhat higher than 
it is derived by other investigators. 

To the best of our knowledge, nobody has investigated the COBE noise 
(A-B) maps in detail. Gorski et al. (1994) used only the mean noise parame- 
ters assumed they are normal. Wright et al. (1994b) used the noise spectrum, 
but the used basis was not complete orthogonal. Bennett et al. (1994) worked 
by the help of a correlation function but only on the part of the sphere. In the 
latter two cases the basis function orthogonality is lost and it is impossible 
to analyze noise spectrum in detail. 

The next result we obtained and that differs from previously published is 
the rejection of the trivial (with n = 1) Harrison-Zel'dovich spectrum on the 
confidence level of 99%. 

As a result there are very unlikely the models with high A-term and a 
more probable are the open models with < 1, if we assume the initial den- 
sity perturbations are determined by the power law index n = 1 (Kamionkowski M., 
Spergel D.N., 1994). Our result supports the existence of a barionic entropy 
models. At the same time there are not ruled out the models having ri > 1 in 
analyzed scales and having more complicated inflation potential (Starobin- 
sky 1992, White et al., 1994 and references therein). 

It must be emphasized that by now the spectrum parameters determining 
accuracy is not enough to draw more deep inferences. Moreover we can 
see the obtained results being very sensitive to the accuracy of the noise 
determining and to the used specific estimation procedure. 

We hope the 4-year COBE investigation circle will be accessible in the 
near future. Unfortunately even the 4-year data will not improve the situa- 
tion drastically. First, it is necessary to increase the instrumental sensitivity 
as a minimum an order. This will be reached in the planned experiment 
" RELICT-2" . So the accuracy of the n estimation will be about 5-7% and will 
be determined the cosmic variance rather than instrumental noise (Sazhin et 
al., 1995). Second, it is necessary to enhance the anisotropy measurement 
angle resolution coupled with a sky coverage increasing (Scott et al., 1994). 
Unfortunately, the modern middle and small scale investigations give us the 
information only about a few point on the sky and so the results are strongly 
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contaminated by the sample variance. 
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Table Caption 



Table |T|. Mean spectral indexes and its variations determined for 
spherical harmonics at different frequency ranges. The first column 

- the frequency ranges where the spectral index a is analyzed. The 
second column - spectral index determined from the most significant 
spectral components with the signal/noise ratio lager than 5, 05. The 
third column - its variation, Acts. The forth column - spectral index 
variance predicted from noise analysis, Aa^. The fifth column - the 
mean weighted spectral index aw for all harmonics with m 7^ for 
/ = 2k, and m 7^ 1 for / = 2k + l. The sixth column - its variance Aaw 

Table |2|. Comparative result of an influence of spatial filtering (i.e. 
excluding the components with m = for I = 2k and m = 1 for / = 2k 
from / = 2 to Z = 25 inclusively) to a signal amplitude for different 
frequencies. The last row shows the signal amplitude after the spatial 

- frequency filtering assuming the Galactic radiation spectral index as 
corresponded to ionized hydrogen. <J(a+b)/2 ^ amplitude, calculated 
as a half-sum of two maps (i.e. signal plus noise), <J(^a-b)/2 ^ noise, 
calculated as a half-difference of two maps, asky ^ signal amplitude 
estimation on the sky. 

Table ^. Spectrum analysis after spatial-frequency Galactic radiation 
filtering. The influence of quadrupole excluding is shown. In addition 
it is shown the results after noise anomalous harmonics are excluded. 
The first row is corresponded to spectrum from / = 2 to / = 22 with 
the abnormal harmonics excluded. The second row - the same but 
in addition the quadrupole is excluded. The third row shows the re- 
sults if all harmonics from / = 2 to / = 25 are used. The forth row 
shows in this case the influence of quadrupole exclusion. cr\_^Q - total 
measured power within given spatial window, cr^_^ - noise power esti- 
mation obtained within the same window from a difference map, (Jsky 

- cosmological signal power estimation, {dgi^y = (t\j^q — a\_^ ), Q2, 
n - estimations of parameters for spectrum (|I]), M - total number of 
analyzed harmonics. In all cases we assume that m 7^ for / = 2k and 
m 7^ 1 for Z = 2A; + 1. 

Table ^. Results of CMB anisotropy spatial spectrum analysis ob- 
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tained by different authors basing on COBE 2-years data. For the 
comparison our result is shown. It is shown how a quadrupole exclud- 
ing does influence on the end result. 
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Table 1: 



frequency 
(GHz) 


"5 


A«5 




aw 




19.2 - 31.5 
31.5 - 53.0 


-2.13 
-2.27 


0.36 
0.28 


0.38 
0.36 


-2.12 
-2.19 


0.69 
0.76 



Table 2: 



frequency 
(GHz) 


spatial 
filtering 


Cr{A+B)/2 


Cr{A-B)/2 


(7 Sky 


31.5 


yes 


641. 


95.92 


634.70 


no 


254.79 


94.52 


239.83 


53. 


no 


191.64 


33.04 


188.80 


yes 


83.32 


32.20 


76.86 


frequency 
" cleaning" 
q; = -2.15 


yes 


50.78 


45.44 


22.67 



Table 3: 



NN 


Quadrupole 
analysis 


2 

'^A+B 


'^A-B 


^Sky 


Q2 
fiK 


n 


M 


1 


yes 


(70.28)^ 


(56.76)2 


(41.44)2 ^ (1458)2 


15.22 ±2.9 


1.84 ±0.29 


446 


2 


no 


(68.92)2 


(56.47)2 


(39.51)2 ± (14.65)2 


15.55 ±3.8 


1.81 ±0.37 


442 


3 


yes 


(82.91)2 


(74.19)2 


(37.03)2 ^ (17.48)2 


18.03 


1.31 


649 


4 


no 








20.3 


1.12 


645 
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Table 4: 



Author 


n 


Qrms—PS 


Quadrupole included 


Bennett et al.(1994) 


-, 40+0.49 


12.811:1 


Gorski et al.(1994) 


1 99+0.43 
-■-•^^-0.52 


i7.oiy 


Wright et al.( 1994b) 


1 Qq+0.34 
^^■'^^-0.39 




this work 


1.84 ±0.29 


15.22 ±2.9 


Quadrupole excluded 


Bennett et al.(1994) 


1 11+0.60 




Gorski et al.(1994) 


1 00+0.53 
-'-•U^_0.59 




Wright et al.( 1994a), Z = 3 - 30 


1 9 c +0.40 
-•-•^^-0.45 




Wright et al.(1994b), Z = 3 - 19 


-'-■^O_o.44 




this work 


1.81 ±0.37 


15.55 ±3.8 
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